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An estimate of the local ISW signal, and its impact on 
CMB anomalies 
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ABSTRACT 

We estimate the local density field in redshift shells to a maximum redshift of z = 0.3, 
using photometric redshifts for the 2MASS galaxy catalogue, matched to optical data 
from the SuperCOSMOS galaxy catalogue. This density-field map is used to predict 
the Integrated Sachs- Wolfe (ISW) CMB anisotropics that originate within the volume 
at z < 0.3. We investigate the impact of this estimated ISW foreground signal on large- 
scale anomalies in the WMAP CMB data. We find that removal of the foreground 
ISW signal from WMAP data reduces the significance of a number of reported large- 
scale anomalies in the CMB, including the low quadrupole power and the apparent 
alignment between the CMB quadrupole and octopole. 
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1 INTRODUCTION 

The possible presence of non- Gaussian features in the low- 
order multipoles of the Cosmic Microwave Background has 
been a recurring theme in cosmology for over a decade. 
Large-scale anomalies of varying significance have been 
claimed, includ i ng th e low power in the C MB quadrupole 
([Hinshaw et al.1 Il996l : ISoergel et all l2003h : the planarity 
of the octo pole and its apparent alignment with the 
quadrupole (Tegmark et al ] (20031 : |de Oliveira- Costa et al.l 
l2004r ); the nort h-south power asym metry about the eclip- 
tic plane (e.g. lEriksen et al.l |2004 ) and the extension of 
the I — 2 and I — 3 alignments t o higher multipoles 
(ISchwarz etHl 120041: ICopi et alJl2006h . defining the Axis 
of Evil' ([Land & Magueiiol I2005T ). A closely related result 
i s the lack of large - angle correlation s in the CMB sky 
dHinshaw et al. 1996; ICopi et al. On smaller scales, 

IVielva et al.l (|2004T ) suggest that the CMB contains a non- 
Gaussian 'cold spot' in the southern hemisphere. There are 
four possible explanations for these CMB anomalies: (i) a 
problem with the data; (ii) an unlucky fluke; (iii) new early- 
universe physics; or (iv) they are due to a foreground signal. 
Here we explore the last of these possibilities. 

The ISW effect is a secondary anisotropy in the CMB 
caused by the passage of CMB photons through evolving 
gravitational potential wells: 
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where to and t LS denote the times today and at last scat- 
tering respectively; x is the position along the line of sight 
of th e photon at time t and <E> is t he gravitational poten- 
tial (Ma rtinez- Gonzalez et al.l fl990). On large scales in a 
Dark Energy dominated universe, CMB photons gain en- 
ergy when they pass through the decaying potential wells 
associated with overdensities and lose energy on passing 
through underdensities. This leads to the large-scale posi- 
tive cross-correlation between large scale structure and the 
CMB that is often used to attem pt detections of the ISW 
effect ( Crittenden fe Turold[l996h . 

The ISW effect is the only significant secondary CMB 
anisotropy at low i. The only other non-negligible signal, due 
to the thermal SZ effect, is approximately an order of mag- 
nitude smaller than the expected ISW signal. Our aim here 
is to predict the local ISW signal to a maximum redshift of 
z — 0.3 and consider the effect of removing this part of the 
foreground signal on the large-scale CMB anomalies. The 
data and method used to predict the local ISW signal are 
described in Section [2] and the effect of this foreground sig- 
nal on large-scale CMB anomalies in explored in Section [3] 
Finally, our conclusions are presented in Section 2] A cosmo- 
logical model with Q m 0.3, Q v 0.7, Q b 0.05, h 0.7 
and erg = 0.75 is assumed unless otherwise stated. 



2 PREDICTION OF THE LOCAL ISW SIGNAL 

In order to predict the local ISW signal, it is first necessary 
to estimate the density field within the volume of interest. 
We make a 2D Wiener reconstruction of the projected den- 
sity field in redshift shells, using photometric redshift data 
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for the 2MASS galaxy catalogue. An account of the data and 
method used for this reconstruction is given in Section |2J"1 
and the estimation of the local ISW signal from this density 
field is described in Section [2?2l 



2.1 Recovery of the local density field 

2.1.1 Data 

The 2MASS dataset used in this paper is described briefly 
below. This material is common to a companion paper about 
attempted detection of th e ISW effect via cross-correlation 
(Fr ancis fe Peacockll2009h and fuller details are given there. 

2M ASS is an a ll- sky survey in the J, H and K s 
bands ([jarrettl 1 20041 ). This near-infrared selection means 
that 2MASS is sensitive to old stars and hence to the most 
massive structures, which makes it well-suited for construct- 
ing maps of galaxy mass density. The final extended source 
catalogue (XSC) contains over 1.6 million objects, over 98% 
of which are galaxies. Photometric redshifts for the 2MASS 
XSC have been generated by matching the 2MASS data 
with optical catalogues from SuperCOSMOS scans of ma- 
jor photographic sky surveys (UKST in the south; POSS2 
in the north) to create a 5-band BRJHK dataset. Details 
of the SuperCOSMOS cata logue construction process are 
given m iHamblv et all (|200lh and details of the photometric 
calibration and phot o- z determination method are given in 
IPeacock et al.l (|2010h . The final photometric redshift dataset 
has an overall rms in z p hot — z of 0.033. 

A magnitude cut K s < 13.8 is imposed on the 
extinction-corrected data to ensure uniformity across the re- 
gion outside the Galactic plane. T he Galactic plane ma sk is 
constructed using the dust maps of lSchlegel et al.l (| 19981 ) and 
excludes all regions with if-band reddening Ak > 0.05. This 
leaves approximately 67% of the sky unmasked. Wh i lst the 
choice of Ak > 0.05 is conservative, lAfshordi et al.l ([2004) 
note that the number density within pixels of faint galax- 
ies (13.5 < K20 < 14.0) drops off as a function of Ak for 
Ak > 0.065. A less conservative mask, for example with 
Ar > 0.1, would still only leave 79% of the sky unmasked - 
at the price of possible incompleteness issues in the analysis. 



2.1.2 Density field reconstruction 

It is possible to attempt a full 3D reconstruction of the lo- 
cal density field with this photometric redshift data set, and 
then use this to estimate the local ISW signal (see I Francis! 
2008). However, this recovered density field is subject to sig- 
nificant radial smearing due to the uncertainty in the radial 
position of each galaxy (approximately 100/i _1 Mpc). The 
extent of this smearing justifies a simpler and more trans- 
parent approach: to recover the projected density field for 3 
thick redshift slices of width Az = 0.1 to z ma x = 0.3, and 
estimate the ISW signal due to each slice. 

The galaxy distribution in each photometric redshift 
shell is known outside the masked region. We use a Wiener 
filtering technique to recover the galaxy density field in each 
shell. When the Wiener filter is applied to all-sky data, it 
recovers the optimal density field in the least squares sense; 
a thorough discussi on of Wiener- fi l ter de nsity reconstruction 
is presented in e.g. IZaroubi et al ] (|l995h . 

Because our dataset is masked at low Galactic latitudes, 



the simple all-sky Wiener approach must be modified. We 
proceed by initially filling the masked region in each of the 
redshift shells with a random Poisson sampling of galaxies 
according to the average number density observed outside 
the plane. This galaxy distribution is then used to deduce 
the projected galaxy density field for the redshift shell. The 
Wiener filter appropriate for each shell is then constructed 
as usual: 



$wf(€) 
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where C gg (£) is the angular power spectrum of galaxies in 
that redshift shell and C no ise(£) describes shot noise. The 
Wiener filter is applied in the first instance to the aim co- 
efficients computed for the expansion of this initial full-sky 
galaxy density field into a basis of spherical harmonics. The 
resultant density field within the masked region is Poisson 
sampled for galaxies once more to yield a revised galaxy 
distribution composed of a sample from the Wiener recon- 
struction within the plane, and the original 2MASS galaxies 
in the unmasked regions. In this way, large-scale features 
in surface density can be extrapolated across the plane in a 
consistent way. The process can be iterated, but this did not 
yield a significant improvement in the final results. We have 
carried out tests of this method on simulated data, and find 
satisfactory recovery of the estimated ISW signal (Francis 
2008). Fig. [T] shows the recovered projected density field for 
each of the three shells. We use HEALPix software (Gorski 
et al. 2005; see also http://healpix.jpl.nasa.gov) to gen- 
erate these maps and to perform harmonic analysis of them. 

2.2 Prediction of the local ISW signal 

Given a projected galaxy density field in a thick redshift 
shell, it is possible to estimate the CMB temperature fluc- 
tuations, AT, due to the ISW effect. Poisson's equation tells 
us the time variation of the potential field in terms of the 
density growth factor in the linear regime, g{p): 



$(a) = ^$( 
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(3) 



Hence, using Poisson's equation, we can integrate through ; 
shell of thickness Ar: 



AT* 



9(a) 



dt 

£(£ + l)c 3 



2 <S>£m(a) dr 



dg ,_v 

-jz( a ) 



9(a) 
9(a) 



9{a y 



I 



®£m dr 



9{a) 



r 2 (a)H(a) S -^Ar, 



(4) 



where S is the projected galaxy density field in the redshift 
shell under consideration and a — (1 + z)~ x with z the red- 
shift at the midpoint of the shell. The bias in each redshift 
shell between the galaxy and matter density fields is approx- 
imated using a linear bias relation S g = bS m and is assumed 
to be independent of scale and redshift in each shell. To de- 
duce the appropriate value of b in each redshift shell, a max- 
imum likelihood fit of the predicted galaxy angular power 
spectrum to the measured galaxy an g ular p ower spectrum 
is performed as in lFrancis &; Peacockl (|2009f ). 

Fig. [1] shows the results of estimating the local ISW ef- 
fect from the Wiener filtered density field in redshift shells. 
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Figure 1. (Left) The 2D reconstruction of the local density field described in Section [2?T1 in three photometric redshift shells: 0.0 < z < 0.1 
(top), 0.1 < z < 0.2 (middle) and 0.2 < z < 0.3 (bottom). The plots show overdensity 5 on a scale -0.6 ^ 5 < 0.6. (Right) The 
corresponding ISW signal in mK computed from the reconstructed density field using equation ([4)>. 



These ISW maps are affected by overlap effects due to the 
uncertainty in the radial position of each galaxy: galaxies 
associated with an overdensity near the boundary of a shell 
will to some extent be spread across both the 'true' shell 
and the adjacent redshift shell. Fig. [2] shows the angular 
power spectrum of the ISW signal in each slice, together 
with the predicted linear ISW signal. For these predictions, 
we have accounted for the smearing effects of the photomet- 
ric redshift data by assuming that photo- z's apply Gaussian 
smoothing along the radial axis with a r = 90/i -1 Mpc. Fig. 
[2] also shows the effect of this radial smearing on the to- 
tal predicted ISW signal to z ma x = 0.3. Overall, the agree- 
ment between observed and predicted power spectra is good, 
which illustrates that the effect of photo- z imperfection is 
relatively minor for £ < 20. 

We note that the estimated ISW signal in the lower two 
redshift slices is slightly larger than expected at £ > 20; this 
is particularly noticeable in the z < 0.1 slice. The nonlin- 



ear Rees-Sciama effect is expected to increas e the power at 
high I, doubling the power at £ ~ 200 fe.g. lCoorav h Shethl 
120021 ). Since the mean depth of the lowest-redshift shell is 
of order one tenth of the distances that dominate the to- 
tal ISW effect, it is plausible that we are seeing some lo- 
cal Rees-Sciama effect. In any case, the main focus of the 
present paper is at larger angular scales. We also note that 
the estimated ISW power is larger than expected for £ < 3 
in the 0.1 < z < 0.2 shell and for £ < 5 in the 0.2 < z < 0.3 
shell, which leads to a larger than average total signal for 
Zmax = 0.3 on such scales. Such small-^ modes are po- 
tentially sensitive to the corrections for mask incomplete- 
ness near the plane, but tests of our method on mock data 
show no tendency to bias the power high in this regime (al- 
though it does increase the already considerable cosmic vari- 
ance). Any discrepancy is relative to our standard model 
of Q m = 0.3, whic h is perhaps on the high-density side 
given current data ([Komatsu et al.l l2009h . Decreasing Q m 
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Figure 2. {Left) The power spectrum of the ISW signal estimated in each of the three redshift shells (dashed lines): 0.0 < z < 0.1 
(bottom), 0.1 < z < 0.2 (middle) and 0.2 < z < 0.3 (top) together with the predicted signal for f2 m = 0.3 (solid lines), allowing for radial 
smearing. (Right) The predicted total ISW power for z max — > oo (upper red solid line) and for z max = 0.3 (middle purple dashed line) 
and that expected for z m ax = 0.3 allowing for radial smearing (lower purple solid line). The green dashed line shows the power spectrum 
measured for the estimated z max = 0.3 ISW signal. 



to 0.25 with other parameters held fixed raises the expected 
quadrupole power for z max = 0.3 by a factor 1.5; our adopted 
normalization of as — 0.75 might also be increased, so that 
the predicted \ow-£ ISW power could plausibly increase in 
amplitude by a factor 2. We will thus provisionally adopt 
the large-scale ISW map as estimated, but will also explore 
the consequences of scaling its amplitude down to the lower 
values expected in our Q m = 0.3 model. 

Fig. [2] also shows the total ISW signal for 
which we can use to estimate the fraction of the total ISW 
signal that we expect to include if z max = 0.3. The square 
root of the ratio of the predicted ISW power for 
to that for z max = 0.3 shows that we can expect to include 
~ 40% of the total rms ISW temperature signal at low mul- 
tipoles, falling to ~ 10% at £ = 30. This variation with 
angular scale is mainly due to the greater importance of the 
high redshift ISW signal at large £. 



3 IMPLICATIONS FOR CMB ANOMALIES 

Having made a prediction of the local ISW signal, we now 
investigate the effect of this foreground signal on a number 
of the large scale anomalies that have been reported in the 
CMB. Many of these CMB anomalies relate solely to the 
£ = 2 and £ = 3 multipoles and, as noted in Section 12.21 
the estimated local ISW power is greater than the ensemble 
average on these scales. We therefore consider two versions 
of the estimated local ISW signal in this analysis of the 
CMB anomalies: the ISW prediction described in Section [2] 
and a rescaled version of this signal chosen to maintain the 
phases of the original signal but to have the power spectrum 
expected of the ISW signal out to z = 0.3 for our standard 
choice of Q m = 0.3. Consideration of this latter ISW signal 
allows us to assess the possibility that any effect on the CMB 
anomalies is solely due to the large low-multipole power of 
the estimated ISW signal. 



We now consider the effect of removing the local ISW 
signal on each of the large scale CMB anomalies. We use the 
year-3 WMA P ILC map (hereafter ILC) and the C MB map 
computed bv lde Oliveira- Costa &; Tegmarkl (|2006r ) from the 
year-3 WMAP data (hereafter OT). We have removed the 
kinetic Doppler dipole and quadrupole contributions due to 
the motion of the Sun with respect to the CMB frame. 

3.1 Low quadrupole power 

Figure [3] shows the £ = 2 and £ = 3 multipoles of the es- 
timated local ISW signal, the WMAP ILC map and the 
WMAP ILC map after removal of the estimated ISW sig- 
nal. We immediately note a significant increase in the tem- 
perature range of the quadrupole map after ISW removal 
and an alteration in the positions of the temperature ex- 
trema for this multipole; in contrast, the octopole is largely 
unchanged. We compute the power of the quadrupole for 
these maps, together with the fore ground- cleaned map of 
Ide Oliveira-Costa fe Tegmarkl (|2006T ) before and after ISW 
subtraction. The results and their significance given the 
WMAP best-fitting model prediction are given in Table [TJ 
We also repeat the exercise for the CMB maps with the 
rescaled version of the local ISW prediction removed. The 
results in Table [1] show that the low value of the quadrupole 
power ceases to be significant when we consider the maps 
with either the original or the rescaled versions of the es- 
timated local ISW signal removed. This is one of the main 
results of this paper, and its broader significance is discussed 
in the concluding section. 

3.2 Quadrupole/octopole alignment 

In order to quantify the alignment between the quadrupole 
and octopole, we have used the maximum an gular disper- 
sion technique of Ide Oliveira-Costa et al.l (| 20041 ) to identify 
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Figure 3. (Top) The quadrupole (left) and octopole (right) in mK of the estimated local ISW effect. (Centre) The quadrupole (left) 
and octopole (right) in mK of the WMAP ILC map. (Bottom) The quadrupole (left) and octopole (right) in mK of the ILC map after 
subtracting the local ISW contribution. 





Map 


Quadrupole Power / /xK 2 


Probability 




ILC 


250.6 


3.5% 




OT 


214.3 


2.4% 




ILC - ISW 


600.7 


19.8% 




OT - ISW 


608.2 


20.1% 


ILC 


- ISW (rescaled) 


372.2 


8.0% 


OT 


- ISW (rescaled) 


358.8 


7.5% 



Table 1. The quadrupole power for the WMAP ILC and OT 
maps before and after subtraction of the estimated local ISW 
signal for z < 0.3, together with its significance given the best- 
fitting year-3 WMAP prediction. This is a one-tailed probability 
of having a quadrupole as low or lower than observed, given the 
distribution of quadrupole power (x 2 with 5 d.f.). We note that 
the quadrupole power is no longer unusually low after subtraction 
of the ISW signal. 



preferred axes for the quadrupole, n2, and octopole, 113. 
Thus we search over directions n to maximise the quantity 



^2 rn m \a£m\ in a coordinate system with n as the direction 
of the z-axis. Our results are given in Table [2] We see that 
there is no significant alignment remaining in either the ILC 
or OT map after removal of either the estimated local ISW 
signal or the rescaled version. 



3.3 Planarity of the octopole 

Fig. [3] shows that the octopole is not particularly affected by 
the subtraction of the estimated local ISW effect. Evaluating 
the statistic of Ide Oliveira-Costa et al. 



a 3 



-3| 2 + 



«33 



Em \ a ^ 



(5) 



as a measure of planarity, we find the results given in Table 
[3] The significance of the '£' statistic is evaluated by con- 
sidering 10000 simulated sets of aim data and determining 
the frequency of such values occurring by chance (see 
de Oliveira-Costa et al. 2004 for further details). The ISW 
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Figure 4. The £ = 2 + 3 multipoles of the year-3 ILC map (left) and the ILC map after subtraction of the estimated local ISW signal 
(right) in mK; the black line shows the ecliptic plane. We note that the power asymmetry persists, with the largest temperature extrema 
seen south of the ecliptic plane, but that the ecliptic no longer lies along a node line of the temperature distribution. 





Map 


n 2 • n 3 


Separation 


Probability 




ILC 


0.9991 


2.4° 


0.09% 




OT 


0.9881 


8.9° 


1.2% 




ILC - ISW 


0.7548 


41.0° 


24.5% 




OT - ISW 


0.6712 


47.8° 


32.9% 


ILC 


- ISW (rescaled) 


0.8888 


27.3° 


11.1% 


OT 


- ISW (rescaled) 


0.8183 


35.1° 


19.2% 



Table 2. The scalar product of the preferred axes of the 
quadrupole and octopole together with their angular separation 
and the one-tailed probability of such a good separation occur- 
ring by random alignment. Results are shown for the ILC and OT 
maps before and after removal of both the estimated local ISW 
signal and the rescaled version. 





Map 


't' value 


Probability 




ILC 


0.9213 


17% 




OT 


0.9269 


15% 




ILC - ISW 


0.9841 


1.6% 




OT - ISW 


0.9653 


5.0% 


ILC 


- ISW (rescaled) 


0.9661 


4.9% 


OT 


- ISW (rescaled) 


0.9524 


8.0% 



Table 3. The t value for the octopoles of each of our usual maps 
together with the one-tailed probability of such an extreme value 
occurring by chance, as determined from 10000 Monte-Carlo sim- 
ulations of a3m coefficients. Removal of the estimated local ISW 
effect actually makes the octopole more planar. 



subtracted maps have octopoles that are slightly more pla- 
nar than the corresponding ILC or OT maps, i.e. removal of 
the estimated local ISW effect seems to make the octopole 
more planar. The probability of such planarity occurring by 
chance varies depending on the CMB map which is used, 
with probabilities of ~ 2 — 5% for the full local ISW sub- 
traction and 5 — 8% for the rescaled ISW subtraction. 



3.4 North-South asymmetry 

Fig. [4] shows the sum of the quadrupole and octopole for the 
ILC map, before and after ISW subtraction. An asymmetry 
in power between the hemispheres north and south of the 
ecliptic remains after removal of the estimated ISW signal, 



but the ecliptic no longer lies along a null of the tempera- 
ture distribution. The significance of the power asymmetry 
and the coincidence of the ecliptic with a node line across 
1/3 of the sky was found bv lSchwarz et all (|2004l ) to be un- 
usual at the 95% level. We expect that the breaking of this 
coincidence with the ecliptic that occurs after ISW removal 
would reduce the significance of the anomaly much further. 



3.5 Cold spot 

The 'cold spot' is centred at £ = 207.8°, b = —56.3° and 
has sc ale ~ 10°. This feature was identified bv IVielva et al.1 
using Spherical Mexican Hat Wavelets in the year- 
1 W MAP CMB data a nd confirmed in the y ear-3 WMAP 
data (jCruz et al.ll2007h . iRudnick et al.1 (|20Q7h identified the 
region of the cold spot with an underdensity of NVSS radio 
sources and therefore suggested that the cold spot is due 
to the IS W contribution of a lar ge void at redshift z ~ 1. 
However, [Smith &; Hutererl (|2Q10h find no evidence for such 
an underdensity of NVSS sources. We have evaluated the 
CMB temperature in the vicinity of the cold spot before 
and after the subtraction of the estimated local ISW signal; 
the pixels in the vicinity of the cold spot become on average 
7/xK warmer after ISW removal. This changes the formal 
significance level of the feature from 1.85% to 2.5%. 

If the cold spot were due to the ISW contribution of a 
large void at z ~ 1, removal of the estimated local ISW 
signal slightly reduces the size of the necessary void to 
around 95 % of the current est imate of about 105/i _1 Mpc 
at z ~ 1 (jRudnick et al.ll2007h . A void of this size would 
remain anomalously large. 



4 CONCLUSIONS 

We have estimated the local ISW signal to z = 0.3 based 
on large scale structure data from 2MASS and SuperCOS- 
MOS. On subtracting this predicted foreground from CMB 
data, we find that many large-scale CMB anomalies are al- 
leviated. In particular the quadrupole power is raised and 
the alignment between the £ = 2 and £ = 3 multipoles is 
broken. We have checked that these results are not sensi- 
tive to the precise amplitude of the low multipoles of the 
estimated ISW signal by considering in addition a rescaled 
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version of the ISW map, chosen to have the same phases as 
the original estimation but a power spectrum equal to that 
expected for the local ISW effect with Q m = 0.3. On sub- 
traction of the estimated ISW foreground, we also find that 
the co-incidental alignment of the ecliptic with a node line 
of the £ = 2 + 3 multipoles of the CMB is broken (although 
the north-south asymmetry in the signal remains) and that 
the temperature in the vicinity of the non-Gaussian cold 
spot is slightly raised (although the change in significance 
is not substantial). Of the large-scale anomalies considered, 
the planarity of the octopole remains most unusual, and in 
fact becomes more significant on removal of the ISW signal 
- but ~ 5% significance is hardly compelling evidence for 
any anomaly. 

These results thus demonstrate that removal of the lo- 
cal ISW effect can mitigate the anomalous signals that have 
been claimed on large scales in the CMB. Our estimate of 
the CMB signal that originates at z > 0.3 has a character 
that is a better match to the expected Gaussian random 
field. The harmonic-space anomalies investigated here are 
closely related to the issue of lack of large-scale angular cor- 
relations, and one might e xpect that ISW subtrac tion would 
also remove this problem; lEfstathiou et al.l (|2009l ) have used 
the data from the present paper to show that this is indeed 
the case. 

It should be borne in mind that we have only considered 
the ISW signal for < z < 0.3 in this work, which we esti- 
mate to contain around 40% of the total rms ISW signal at 
low multipoles, and therefore there will be further substan- 
tial influence on the observed CMB from the ISW effect at 
higher redshifts: an upper limit of z < 1.3 would yield > 90% 
of the signal. Future imaging surveys will provide photo- 
metric redshifts to this depth (approximately an AB mag- 
nitude limit of r = 24), but extensive s ky coverage will be 
hard to obtain; the Dark Energy Survey (|DES Collaboration] 
l2005h will cover 5000 deg 2 by about 2015, but the full extra- 
galactic sky will p robably have to await results from LSST 
Ivezic et al. l l2008h in perha ps 2020, or ultimately ESA's Eu- 



clid mission (|Laureiid l2009). which is currently scheduled to 
finish in about 2023. 

Our work nevertheless demonstrates that it is important 
to estimate the contribution of ISW secondary anisotropics 
to the CMB before concluding that the large scale in- 
trinsic CMB signal is anomalous or otherwise. In this re- 
spect, we agree with previous workers who have argued 
that foregrounds from large-scale structure might affect the 
CMB anomalies: iRakic et al.l (|2006|) lo o ked a t the nonlinear 
Rees-Sciama effect and Abramo et al.l (|2006l ) the Sunyaev- 
Zeldovich effect. However, it is well established that the ISW 
signal is much the largest of these fo regrounds on the rele- 
vant scales fe.g. lCoorav &; Shethll20oj ). and so ISW removal 
should dominate any change in the character of the CMB 
sky. 

Of course, it can be argued that we have replaced one 
puzzle by another: the intrinsic CMB should be statisti- 
cally isotropic and Gaussian, and the large-angle ISW signal 
should very nearly have the same character. At the low- 
est multipoles, the prediction of the standard cosmology 
is that the power from each of these co mponents should 
be ve ry nearly equal (see Figure 62 from ICoorav &; Shethl 
2002). How could the sum of two uncorrelated Gaussian 
fields reduce the total quadrupole power, and introduce non- 



Gaussian correlations at higher multipoles? Within the stan- 
dard model, such an outcome woul d have to be a flu ke, and 
it is true that some authors (e.g. lEfstathioul 120041 ) claim 
that the statistical significance of the claimed anomalies is 
not overwhelming. 

If we reject the idea of a fluke, then either the intrinsic 
CMB or the foreground would have to be anomalous, unless 
a new mechanism can be found to correlate the se two compo- 
nents. Regarding non-standard foregrounds, iGordon et al.l 
( 2005) showed that it is hard to produce the complete set of 
claimed anomalies by adding uncorrelated signals even if the 
foreground component is made explicitly anisotropic; rather, 
they claim that one would need a multiplicative anisotropy 
in which the temperature variance has a large-scale varia- 
tion over the sky. At higher multipoles than considered here, 
there is accumulating evidence that just such an effect ex- 
ists, with both a dipole and quadrupole modulation (e.g. 
lHanson fe Leldsl 120091 : iHoftuft et aill2009h . If such effects 
are not low-level instrumental systematics, then a radical 
new ingredient in cosmological structure formation is re- 
quired; but in either case, there is no evidence that these 
modulations have any relation to the claimed large-angle 
anomalies. Our position on these is that exotic mechanisms 
for explaining CMB anomalies as a real physical effect (e.g. 
iLuminet etaU l2003h . tend to apply only to the intrinsic 
anisotropics that originate at z ~ 1100, so that we should 
try to isolate these from foreground emission as cleanly as 
possible. We have tried to perform this task in the present 
paper by removing part of the ISW foreground, with the 
result that the estimate of the intrinsic CMB no longer ap- 
pears unusual in its large-scale properties. This reduces the 
motivation for revisions to the standard theory for the ori- 
gin of CMB fluctuations, even if it leaves the door open for 
a more elaborate model that could introduce a correlation 
between the intrinsic CMB and foreground anisotropics. 
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